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1.  Introduction 


In  2003,  Wang  et  al.1  derived  expressions  for  the  Differential  Group  Delay  (DGD)  of  a  randomly 
birefringent  fiber  in  the  Fixed  Modulus  Model  (FMM).  The  DGD  has  both  modulus  and  phase. 
The  FMM  assumes  that  the  modulus  of  the  birefringence  vector  is  a  random  variable.  Wang  et  al. 
presented  analytical  results  with  the  following  assumptions: 

•  The  spin  function  is  periodic. 

•  The  periodicity  length  (p)  of  the  fiber  is  much  smaller  that  the  fiber  correlation  length  (Lf) 
or  p«  LF. 

Later  they  also  generalized  the  FMM  and  presented  the  Random  Modulus  Model  (RMM),  which 
also  includes  the  randomness  in  the  direction  of  the  birefringence  vector.  Now  the  RMM 
equations  could  only  be  solved  numerically. 

In  the  present  work,  the  full  implications  of  the  FMM  have  been  explored  under  the  following 
conditions: 

•  The  p«  Lf  approximation  has  been  relaxed. 

•  A  nonzero  twist  has  been  included. 

•  A  constant  spin  rate  has  been  added. 

We  give  the  analytical  solutions  of  the  exact  FMM  equations  and  present  some  numerical  results 
showing  the  effect  of  different  physical  conditions. 


2.  Theoretical  Analysis 

2.1  The  Model  With  Periodic  Spin  Function 

The  starting  point  is  the  well-known  vector  equation  describing  the  change  in  the  Jones  local 
electric  field  vector  A  (co,  z)  with  the  angular  frequency  co  and  distance  z  along  a  twisted  fiber. 


dAx(z)l  dz 

=  -(A  P) 

dA2(z)/dz 

2 

0 

2/00) 


e2/0(z) 


0 


Mz) 

A  (z) 


(i) 


I  Wang,  M.;  Li,  T.;  Jian,  S.  Analytical  Theory  for  Polarization  Mode  Dispersion  of  Spun  and  Twisted  Fiber.  Optics  Express 
2003, 11  (19),  2403-2410. 
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Here  A (3  (co)  is  the  birefringence  and 


oc 

0(z)  =  —  sin(/7z)  (2) 

n 

is  the  periodic  spin  profile  function  with  spin  magnitude  a0  and  angular  frequency  of  spatial 
modulation  r/ . 


The  boundary  conditions  are 


A,  (0)  =  1,  dA{  (0)  /dz  =  0 

(3a) 

A2  (0)  =  0,  dA2  (0)  /  dz  =  i(Ap  /  2) 

(3b) 

Let  s  =  r/z  be  a  dimensionless  variable.  We  use  (d  /  dz)  =  rj(d  /  ds)  to  rewrite  equation  1. 


~AsO)' 

=  ia 

0 

2  tcs  ins 

e 

~AO) 

AW 

— 2/csins 

e 

0 

AA) 

The  subscripts  denote  differentiation  ( Als  =  dAu  / ds,Als  =  dA2s  /  ds ).  Also,  a  =  (A/>/2ij)  and 
c  =  ( a0  / if)  are  dimensionless  constants. 


We  express  all  parameters  in  terms  of  the  lengths  given  as  beat  length  (LB  =  In  I  A/3 ),  spin 
period  ( A  =  2n  /  77 ),  and  coupling  length  ( l0  =  2n  /  a0 ).  Then  we  can  write  a  =  A  /  2 LB  ,c  =  LB/l0. 


The  new  boundary  conditions  are 


A(0)  =  i,A,(0)  =  o 

(5a) 

A2  (0)  =  0,  A2s  (0)  =  ia 

(5b) 

These  equations  (equation  1  or,  equivalently,  equation  4)  do  not  have  analytical  solutions. 

In  the  perturbative  approximation  (see  appendix),  an  analytical  result  has  been  derived 
earlier.1  In  the  present  work,  we  derive  analytic  solutions  by  approximating  the  sine  function  by 
linear  segments  and  compare  them  to  the  perturbative  solutions  for  the  same  segments. 

2.2  Linear  Segment  Approximation  to  the  Periodic  Spin  Function:  Analytical  Solutions 
for  the  Jones  Amplitude  Equations 

2.2.1  The  Model 

The  periods  of  the  straight  line  segments  shown  in  figure  1  approximate  the  periodic  sine 
function.  A  single  period  with  three-segment  approximation  is  shown. 
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0(s) 

nil 


In 

3  nil 

Figure  1 .  The  three-segment  approximation  to 
the  periodic  sine  function. 

The  field  amplitudes  for  a  given  segment  satisfy  the  following  equations. 

A/"V)1  •  r  0  e2W-MJ  A(m)(s)l  (6) 

—  Id 

A2s(m)(5)j  L^2i*mW  o  JLA(m)(^)_ 

The  superscript  and  subscript  m  both  indicate  the  segments  for  which  the  coupled  equations  hold. 
The  limits  of  segments  follow. 

We  require  that  the  endpoints  of  9m(s )  should  be  the  same  as  that  of  the  sine  function  spin  profile 
0(.s)  \spin=  c sin  s  for  all  segments.  Define  c  =  (2c  /  n)  so  that  the  end-point  conditions  for 
segments  hold. 

For  n  =1,  Segment  I:  0<s</r /2, 

9x(s)  =  cs 

Q(s  =  0)\spi  =0  =  ^  =  0), 

G(s  =  7T/ 2)  \spin=  c  =  ex(s=n!T) 

For  n  =2,  Segment  II:  7i 1 2<  s  <3k 1 2, 

92{s)  =  —cs  +  2c 

0(5  =  jcl 2)  \spi=  c  =  92(s  =  n! 2), 

Q(s  =  37r/2)\spi=-c  =  92(s  =  37r/2) 
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For  n  =3,  Segment  III:  3n/2<s<2n, 
8}(s)  =  cs  -4c 

©0  =  3n /  2)  \spi=  -c  =  03(s  =  3a/2), 
0(5  =5k  12)  \spin=  c  =  03(s  =  5n / 2) 


2.2.2  The  General  M-Segment  Solutions 

The  solutions  for  the  mth  segment  have  the  following  general  form: 


me 


■iOm(s)  4  On) 

L1 

(in) 


0) 


(m)  ( m ) 

<2j  +  lb , 


(m)  (m) 

a2  +  jc>2 


fe,/,*: 


2  1  y»,“i 

(m) 


eosr/,,,5 
sin  ry,„5 


(7) 


with 

<lm2  =a2  +0m\s),0mls  =  d0m(s)/ds 


The  exact  solutions  for  the  coupled  equations  in  one  segment  are  related  to  those  in  the  previous 
adjacent  segment  by  the  following  chain-relations  among  the  coefficients. 

Define  u  =  (qm  ]  / qm),v  =  (8m/s  -0m_lls)l qm,  and  then  the  chain-relations  are  given  by 


(m) 


(m) 


(m) 


( m ) 


h 

t3 

0 

o' 

?4 

^2 

0 

0  " 

^2 

h 

0 

0 

+  W 

—  t3 

h 

0 

0 

0 

0 

h 

0 

0 

*4 

0 

0 

^2 

^4  _ 

0 

0 

—  t3 

h 

~  0 

0 

_^2 

-f4l 

0 

0 

h 

F 

+  V 

^2 

^4 

0 

0 

_~t  1 

—  a3 

0 

0 

(w— 1) 

1 

(m- 1) 


(/w— 1) 


(aw— 1) 


Here,  the  matrix  elements  are 

tx  =  cos  q 
t3  =  sin  q 


m- 1  ^ m  I 


C0S  4,  A- 1 >  r2  =  COS  CLn-\Sm-\  Sm  <hnSm- 1 » 

cos  .  *4  =  sin  q,„  ,-V,  sin  cym5m_, . 


(8) 


The  matrix  chain-relations  can  be  written  compactly  by  expressing  the  4x4  matrices  as  outer 
products  (denoted  by  the  symbol  0 )  of  two  2x2  matrices  as 
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h  h 
1 2  h 


h  h 


0  -1 

1  0 


1  0  fl, 


0  l  a. 


2.2.3  The  Specific  3-Segment  Solutions 

The  details  about  solutions  for  3-segments  follow. 
•  Segment  I:  0<s<nl2 


The  equations  are 


\s\s) 

W) 


2i*wl  A?\s) 

0  A2(1)(s) 


,-2  <640) 


The  boundary  conditions  are 


[A®  (^  =  0)  J = 1,  [A, (1)  (s  =  0)  J = 0 

[^(s  =  0)J=0,[a2/1)(^  =  0)\=ia 


(Ha) 


(Hb) 


=  (c/q)=  1  + 


then  the  analytical  solutions  are  similar  to  those  given  in  section  2.2 


e-'SsA^(s) 


(q/a)eicsA2u(s) J  [0  i  J|_sin#s 


1  -  in~f  cos  qs 


Comparison  with  general  expression  gives  the  following  coefficients: 

a*l>  =  1,  b["  =  0,  a2(1)  =  0,  b 2(1)  =-n  ( 1 

For  calculating  Polarization  Mode  Dispersion  (PMD)  Correction  Factor  (PCF),  the  amplitudes 
have  to  be  differentiated  with  respect  to  to,  which  will  be  denoted  by  subscript  to.  Some  useful 
relations  needed  for  this  are 


i  t  \ 

a  a 

—  —  =  n 
dcoyq  y 


> aco  = 


da  y  d(AB )  ( a  Y  a, 

—  =  — ,  y  =  — — — ,  nm  =  -n  —  — 
dco  2 1)  dco  yqjyq 


’dco=a 
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Then  we  can  write 


(< qla)e~rcsAjl\s ) 
eCcsAj\s) 


\<l  J 


P?+iPi'> 

pf+ipf 


P  3 


U) 


Pi 


(1) 


'  l>4  5 


'/A 


a) 


cos  qs 
sings 


(16) 


where 

P =°,^1)  =-nqs,p?  =-qsP*  =n  (17a) 

Pf  =OjJPe1)  =(1-«2)^jP71>  =0?  Ps)=n2  (17b) 

Some  interesting  relations  are  found  as 

A/?  =  (Ajsq!  A)Vl  —  n2 ,  z  =  (AI2tt)s,  a0  =  (2k2  q2  /  A)nVl  -  n2  (18) 


•  Segment  II:  n!2<s<  3k  1 2 
The  equations  are 


Als(2)(s)l  .  I"  0  e2i02(s)l Af\s) 

—  Id 

_As(2,(s)J  [ e~1Wlis)  0  Jl4(2)(^)_ 

(19) 

The  boundary  conditions  are 

II 

w 

fS 

II 

c7 

"E* 

ii 

w 

(20a) 

II 

cn 

<n 

ii 

II 

w 

^T. 

(20b) 

Similar  expressions  exist  for  A,  (s) .  Using  the  chain-relations  with  n  =  2,  the  analytical 
solutions  are  obtained. 


e-i{-cS+2c)A(2  )(j) 

l-n2  +n2  cos  mi  -  in  sin  mi  n(n  sin  mj  +  i  cos  mi) 

cos  qs 

_{qla)ei(-~cs+2c)A™(s)_ 

-  77(1-  cos  7vq)  nsmmi  +  i 

sings 

(21) 


The  to  -differentiated  amplitudes  are  found  as 


(q/a)e 

i(cs—2c) 


—i(cs—2c) 


(2) 


n® 


(s) 


A®<2)(^) 


e  ; 


a<2> +?>22) 
Ps2>  +iPe2) 


P(32>  +iP42) 
Pi2)+iPs2) 


cos  qs 
sin  qs 


where 

p\ 2)  =  77 2  {2(1  -  cos  mi )  —  mi  sin  mi  +  qs  sin  7iq\ 

=  77 (sin mi -mi cos  mj  +  qs cos  m/\ 
p\  =  77' (-2 sin mi  +  mi cos  nq)  -{l-tr  +n~  cos  mf)qs 
p{2)  =  /7  { —  (cos  mj  +  /Tgsin  mj)  +  g.vsin  mj\ 


(22) 


(23a) 

(23b) 

(23c) 

(23d) 
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pp  =  n{(l  -  2n){\ -  cos  nq)  -  (1  -  n2)nq  sin  nq  +  (1  -  n2  )qs  sin  nq\ 
p(2)  =(1  -n2)qs 

p(2)  =  n\—  (1  -  2n2)  sin  nq  +  (1  -  n2)nq cos  Tvq  +  (1  -  n2)(l  -  cos  mi)qs\ 


pP  = 


(23e) 

(23f) 

(23g) 

(23h) 


Segment  III:  3n/2<s  <2n 


The  equations  are 


A/3)(s) 

W) 


A0)(J) 


0  I  A2  (s) 


The  boundary  conditions  are 


A™  (s  =  3n/2)\  =  [a/3)  (s  =  3tt/2)J 
Als(1)  (s  =  3n  1 2)\  =  [a1s(3)  (^  =  3tt/2)J 


(25a) 

(25b) 


Similar  expressions  exist  for  A<3)(^)-  Using  the  chain-relations  with  n  =  3,  the  analytical 
solutions  are  obtained. 


e-'(“-4c)A<3)U) 

(q  /  a)e'(“_4c)  A(3)  (s) 
fl-n2  +n2  cos  jvq 


n2  sin  2 jvq 

[+  (1  -  n  )(sin3;zg-sin;zg) 
«(cos  7tq  -  cos  3 nq) 


n2  sin  2 7iq 

\n 2  cos  2 nq 
-in< 

[+  (1  -  n2  )(1  +  cos  3 nq  -  cos  nq) 
n{sin  nq  - sin3nq) 


+  in 2  (sin  3  nq  -  sin  2;zg  -  sin  ;zg)_  _+  /{l  -  n2  +  n2  (cos  nq  +  cos  2  nq  -  cos  3  nq) } 


cos  qs 
sings 


The  co  -differentiated  amplitudes  are  found  as 


(q/^e-^Aj^is)  JaPpP+ipP  p? +ip?T  cos  qs 
e^s-4c)Af\s)  {  q  )[pP  +  ipT  PP  +ipP JLsing^. 


where 


p\  =  2n^(l-cos2nq-nqsm2nq)  +  n  gssin2yzg 

-  3 n2  sin  2 nq  -  (1  -  3/r  )(sin  3 nq  -  sin  nq) 

pP  =n  +  nq{2n2  cos  2nq  +  (1- n2 )(3 cos 3 7vq- cos 7iq)} 

-  {i n 2  cos  2 nq  +  (1  -  n2)(l  -  cos  7iq  +  cos  3 nq)}qs 


(28a) 


(28b) 
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✓3-v  2  2  2 

p3  =  -2/r(sin  2m/  -  m]  cos  2nq)  -  (1  —  n+  n  cos  2 jvq)qs 

3 n~  cos  2 to]  +  (1  -  3/r  )(1  -  cos  7iq  +  cos  3 7rq) 

+  7jq\2n2  sin  2 7jq  +  (1  -  n2  )(3sin  3 7jq  -  sin  7iq)} 

-  \n 2  sin  2to]  +  (1  -  n 2  )(sin  3mj  -  sin  7iq)  }qs 

p(53)  =  n(\  -  In2  )(cos  370]  -  cos  70j)  +  n(\-n2  )(sin  3 70]  -  sin  to])70] 
+  n(  1  -  /72)(sin^r  -  sin37cq)qs 


P?  =  n 


(28c) 

(28d) 


(28e) 


P(3)  =  (1  -n1)qs  +  n2 


(2  -  3 n2  )(sin  TO]  +  sin  2 70]  -  sin  3nq ) 

+  (1  -  n2  )(3cos  3  to]  -  2  cos  27iq  -  cos  70])70] 
-  (1  -  n2  )(1  -  cos  70]  -  cos  2 TO]  +  cos  3 TO])qs 


p{3)  =  n(l  -  2779)(sin3^  -  sin^g)  +  77(1  -  772)(cos7zg  -  3  cos  3  to])to] 
+  77(1  -  77 2  )(cos  370]  -  cos  TO])qs 


PI 


(3) 


=  77 


+  77~ 


(2  -  3t72  )(1  -  cos  TO]  -  cos  +  cos  3^g) 

+  (1  -  772)(3sin3^  -  2sin  2to]  -  sin  to])to] 
+  (1  -  77  2  )(sin  to]  +  sin  2 70]  -  sin  3 70j)qs 


(28f) 


(28g) 


(28h) 


2.3  Calculation  of  PCF 

The  sum  of  squares  of  the  co-differentiated  amplitudes  is  similar  to  power  and  can  be  calculated 
by  the  following  expression: 


(m) 


(s) 


+ 


A. 


(m) 


2  co 


(s) 


-(V) 


+  ( To) 


(a0j  1  qy 

a -n2){(pr)2+(Jpf))2+(Jpr)2+(Jpr)2}' 

+(^m>)2 +cJpr))2 +(^7m))2 +(^sm))2 

(1  -  77  2){(/7i('"')2  +{P^)2  -(p)"'1)2  -(PT')2} 
+  (P^)2  +  (p(6m))2  -(P^)2  ~(p(sm))2 


(m)  (m)  .  (m)  (m)  I  ,  (m)  (m)  ,  (m)  (m) 

iyl  F 3  ^  F2  FA  ) _r  ^5  FI  ”r  F6  F 8 


+  [(l-  772){ 

Here,  777  (=1,2,  3)  refers  to  segments  in  sequential  manner. 


cos  2qs 
]sin  2^^ 


(29) 


For  calculating  the  normalized  PCF,  we  need  a  similar  expression  for  unspun  fiber: 


[iA..wr-+iA,„Mi: 

(a(Jq)2 


nspun- fiber  ^^^2 


(30) 
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Then  the  expression  for  the  PCF  becomes 


PCF(m)(s)  = 


Aj%) 

2 

+ 

4,(m)(^) 

2 

!  A^)  i2  + 1 

(s')  1^ 

^  '  Junspun- fiber 

1/2 


(31) 


The  left-hand  side  of  equation  31  is  a  function  of  parameters  n  and  q  and  argument  s  .  In 
general,  the  expressions  are  quite  complicated,  but  for  the  first  segment  the  PCF  is  easily 
calculated  and  is  given  by 


PCFa\s ) 


' 

(  ■  \ 

2" 

1  -n2\ 

1- 

sint/s 

l  Qs  ) 

(32) 


3.  Numerical  Results 


The  physical  constants  ( A/3  ,a0 .// ),  or  equivalently  ( LB ,  l0 ,  A  ),  and  the  parameters  ( n,  q ) 
appearing  in  the  PCF  expressions  are  related  by 


q  = 


1  + 


v4Luy 


,  n  = 


-,-1/ 


1  + 


\4LbJ 


(33) 


We  show  results  for  sets  of  parameters  in  two  limits  to  emphasize  the  difference  between  the 
exact  and  perturbative  calculations. 

3.1  The  Small-Q  Limit  (A  <LB) 

As  shown  in  table  1  and  figures  2  and  3,  in  this  limit,  two  sets  of  parameters  were  chosen  to  get 
small  //-values  (less  than  1).  This  corresponds  to  beat  length  being  larger  than  spin  period. 


Table  1.  PCF  vs.  z  plots  with  small-q-limit  parameters. 


Parameters:  A ,LB,l0 
(m) 

Values 

( n,q ) 

Comments 

(1,12,1) 

(0.9978,  0.6379) 

A  «Lb 

(1,5,1) 

(0.9879,  0.6444) 

A  <Lb 
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PMD  Change  Factor  PMD  Change  Factor 


PCF  vs  s 


Dimensionless  distance  (s) 


PCF,  pert 
- PCF,  ex 


Figure  2.  The  PCF  curve  for  a  perturbative  limit  with  A  =  1  and  LB  =  12.  The  curves  for  exact 
and  perturbative  calculations  are  almost  identical. 


PCF  vs  s 


Dimensionless  distance  (s) 


PCF,  pert 
- PCF,  ex 


Figure  3.  The  PCF  curve  for  a  perturbative  limit  with  A  =  1  and  LB  =5.  The  curves  for  exact 
and  perturbative  calculations  are  almost  identical.  Note  that  after  s  =  5,  the  two  start 
diverging  a  little. 
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3.2  The  Large-Q  Limit  (A >LB) 


As  shown  in  table  2  and  figures  4  and  5,  in  this  limit,  two  sets  of  parameters  were  chosen  to  get 
large  ^-values  (much  larger  than  1).  This  corresponds  to  beat  length  being  smaller  than  spin 
period. 


Table  2.  PCF  vs.  z  plots  with  large-q-limit  parameters. 


Parameters:  A ,LB,l0 
(m) 

Values 
( n,q ) 

Comments 

(5,1,1) 

(0.7864,  4.0475) 

A  >  Lb  (physical  nonperturbative  limit) 

(12,1,1) 

(0.7864,9.7139) 

A  »  Lb  (physical  very  nonperturbative  limit) 

PCF  vs  s 


- PCF, ex 

PCF, pert 


Figure  4.  The  PCF  curve  for  a  nonperturbative  limit  with  A  =  5  and  LB=  1.  The  top  and  bottom 
curves  show  exact  and  perturbative  calculations,  respectively.  Note  that  perturbative 
approximation  underestimates  the  PCF  in  this  regime. 
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PCF  vs  s 


PCF,  exact 
PCF, pert 


Figure  5.  The  PCF  curve  for  a  nonperturbative  limit  with  A  =  12  and  LB  =  1 .  The  top  and  bottom 
curves  show  exact  and  perturbative  calculations,  respectively.  It  is  seen  that  perturbative 
approximation  underestimates  the  PCF  in  this  regime. 


4.  Conclusions 


It  was  shown  through  these  calculations  that  the  perturbative  approximation  made  in  Galtarossa 
et  al.2  has  limited  validity  compared  with  an  exact  calculation.  The  three-segment  approximation 
given  here  can  be  extended  to  any  number  of  segments.  At  the  next  level,  the  segment 
expressions  can  be  derived  for  any  given  profile  of  the  spin  function.  The  exact  analytic 
expressions  allow  a  physical  understanding  of  the  limits  of  the  approximations  employed  earlier. 


“  Galtarossa,  A.;  Palmieri,  L.;  Pizzinat,  A.  Optimized  Spinning  Design  for  Low  PMD  Fibers:  An  Analytical  Approach. 
J.  Lightwave  Technol.  2001, 19  (10),  1502-1512. 
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The  perturbative  approach  is  based  on  the  following  assumptions: 

•  The  coupling  between  the  polarization  states  is  so  small  that  the  equations  become 
decoupled. 

.  The  top  component  is  constant  ( At<m>  =  1, m  =  1,2,3 ),  and  only  the  second  component 
changes. 


•  The  boundary  conditions  remain  unchanged. 


Under  these  assumptions,  the  dimensionless  constant  q  becomes  c  ,  which  is  related  to  the 
physical  lengths  as 


c 


'a' 

\loJ 


(Al) 


The  new  equations  and  their  solutions  take  the  following  form: 

7T 

•  Segment  I:  0  <  s  <  — 


Perturbative  equations: 


Solutions: 


1 

>> 

1 _ 

=  ia 

0 

lies 

e 

T 

g-2 

0 

0 

A2U)(s)  = 


(-) 

v£j 


ie  1CS  sines 


(A2) 


(A3) 


The  sum  of  squares  of  the  to-differentiated  amplitudes: 


r'Aj\s)2  +  Ajl\s) 


2  3 


(a^/c) 


=  ^-(l-cos2c,s')  =  sin2 


cs 


(A4) 


pert 


So 


PCF  (s)  pert  = 


4.%)  +a2J1\s) 


pert 


I  Alft)(s)  h  + 1  A^s)  \~  unspuH_fiber 


1/2 


sines 


cs 


(A5) 
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•  Segment  II:  —<  s  <  — 
2  2 


Perturbative  equations: 


a/2W 

V2)W 


=  la 


2i(-cs+2c ) 


0  c 

-2/(— cs+2c)  q 


(A6) 


Solutions: 


=  e 


i(cs— 2c) 


f-1 

yCj 


\-  (1  -  cos  2c)  cos  c  s  +  (sin  2c  +  /)  sin  c  s] 


(A7) 


The  sum  of  squares  of  the  co-differentiated  amplitudes: 


(2) 


1  CO 


(s) 


+ 


(2) 


A2  co 


(s) 


]_ 

pert  ^ 


\c  ) 


[(3  -  2  cos  2c)  +  (cos  4c  -  2  cos  2c)cos  2c 
j  +  (sin  4c  -  2sin  2c)sin  2c s  ; 


Expression  for  PCF  is  obtained  as  before. 
3tt 

•  Segment  III:  —<  s  <  2n 


Perturbative  equations: 


A/3)(5) 

A,(3)(s) 


=  Id 


2i(cs-4c ) 


0  c 

y-2i(cs-4c)  q 


Solutions: 


A  _  g'(-"+4<;) 


f-1 

v? 


[(-1  +  cos2c-i-cos4c-cos6c)] 

<  (cos  c  s 

[+i(-sin2c-sin4c  +  sin6c) 

f(sin2c  +  sin4c-sin6c)  ] 

+  (  (sines 

[+  /(I  +  cos  2c  +  cos  4c  -  cos  6c)J 


The  sum  of  squares  of  the  co-differentiated  amplitudes: 


aj3)w|2+|a;3)w|21  =^(V 

J  pert  Z\C  J 

The  PCF  can  be  calculated  as  before. 


(5 -4  cos  4c) 

+  (2  cos  10c-  cos  8c  -  2  cos  6c)cos  2c  s 
[+  (2  sin  10c  -  sin  8c  -  2  sin  6c)sin  2c s 


(A8) 


(A9) 


(A10) 


(All) 
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